Skin Friction in Simple Wall - Bounded Shear Flows in the Large Reynolds Number 

Limit. 



Victor YakholE 



O 

(N 

0\ 



s : 

I ■ 
=3 ■ 

O 

• I— I ■ 

Oh: 



> 

oo 
m 

O 

o\ 

O 



X 



Department of Mechanical Engineering, 
Boston University, Boston, MA 02215 
(Dated: today) 

A global approach to analysis of fully developed turbulent flows in pipes/channels and zero pressure 
gradient boundary layers is proposed. A new dynamic definition of the boundary layer thickness 
5{x), where x is the distance to the plate origin, is proposed. The Coles - Fernholtz empirical 
correlation for skin friction A — |^ oc l/ln^(5(a;) and S{x) cc a;/ln^(^) are derived from the 

Navier-Stokes equations in the limit Rex — ^ oo. Here r„ and Uo are the wall shear stress and 
free stream velocity, respectively. The theory is formulated as an expansion in powers of a small 
dimensionless parameter Q in the limit x — oo. 



The law of variation of skin friction with Reynolds num- 
ber in turbulent wall flows is one of the oldest riddles of 
physics of turbulence. In addition to the difficulties as- 
sociated with a general problem of strong isotropic tur- 
bulence, the presence of solid walls is responsible for ap- 
pearance of two different characteristic velocities. The so 
called friction velocity, reflecting properties of the near- 
wall sublayer, is defined as ul = v\ ^^Qy'' \waii, so that for 
the dimensionless distance to the wall y+ = = 0(1), 
the ratio U+ ~ U /u, is independent upon Reynolds num- 
ber. In fully developed pipe/channel flows, the parame- 
ter can be expressed in terms of a prescribed constant 
pressure gradient (or gravity) and the friction factor re- 
lates mean velocity to a driving force. In the vicinity of 
the centerline {y/H w 1), the velocity U{y) « Ucl must 
be found as a solution to dynamic equations of motion. 
In the zero-pressure-gradient boundary layers with exter- 
nally prescribed free stream velocity Uq, this parameter 
is related to the shear stress at the plate and, in addtion, 
to the boundary layer thickness S(x), which in this case, 
depends upon distance to the origin x. It is the interplay 
of these two characteristic velocities which makes theo- 
retical evaluation of velocity profiles U{y) a very difficult 
problem. 



Since the skin friction X — — 



(^)2^ all previous 



calculations heavily relied on a detailed knowledge of the- 
oretically (and experimentally) uncertain function U{y) 
needed for calculation of mean velocity U . The analysis 
of pipe/channel flows is typically based on an assumed 
scaling relation for velocity represented in the "inner" 
and "outer" regions of the flow as U{y) = u^,f{y+) and 
U{y) = UcL — UQg{jj), respectively [1]. The parame- 
ters Uh. and Mo are corresponding characteristic veloci- 
ties. Then, different matching conditions applied in the 
" overlap" region lead to different shapes of velocity pro- 



file U{y). To obtain the functional form of U{y) from 
a systematic local theory, one must derive an expression 
for the distribution of the Reynolds stress Tx,y{y) which 
is equivalent to solution of a proverbial " turbulence prob- 
lem" . Therefore, at the present time, scaling of skin fric- 
tion with Reynolds number remains an unsolved problem. 
In a recent paper, assuming the logarithmic velocity pro- 
file across a zero-pressure-gradient boundary layer, Nagib 
et al [2] developed an asymptotic expansion, leading to 
the so called Coles-Fernholtz relation [3]: 



A cx 



In^ Re 



(1) 
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widely accepted as an accurate large Reynolds number 
asymptotics. While this work is based on a solid mathe- 
matical analysis, its starting point, logarithmic profile, is 
an assumption not following the Navier-Stokes equations. 
In the present paper we present a simple global ap- 
proach, leading to the expression (1) for the skin friction 
and 5{x) oc a;/ln^ ^ for the thickness of zero-pressure- 
gradient boundary layers not relying upon any informa- 
tion about local features of wall flows. The theory is 
based on the following concept: Isotropic and homoge- 
neous turbulence can be characterized by a single di- 
mensionless parameter, called Reynolds number: Re — 
UrmsL/v — T' oo whcrc M^ms = ""^ ^^'^^ ^ IS an integral 
scale at which energy is pumped into the system due to 
external forcing or large-scale instability. Various renor- 
malization procedures based on perturbation expansions 
in powers of this ("bare" ) Reynolds number, led to effec- 
tive or renormalized, viscosity vt ~ UrmsL, widely used 
in engineering turbulence modeling. The main outcome 
of the method is a reformulated perturbation series in 
powers of "dressed" or renormalized Reynolds number 
Rex = UrmsL/vT = 0(1), for which no resummation 
method has been developed. It is the lack of a small pa- 
rameter approaching zero in the limit Re — > oo which 
makes this problem so hard. 
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In this respect, the situation with wall flows is at least 
as difficult and evaluation of the energy spectrum and 
scaling exponents of structure functions is an unsolved 
problem. However, in this case, as Re — > oo, the global 
dimensionless parameters Ut,/U — >■ and ^^^^ — >■ 
are small and can be used for construction of the well- 
behaved perturbation expansion leading to prediction of 
global properties of wall flows. This is the main goal of 
this paper. 

Channel/Pipe flows. First, we consider a steady fully 
developed flow between two infinite plates separated by a 

gap yg ~ 2H, so that H < y < —H. (The ccnterline is at 
y = 0).The flow is driven by the pressure gradient ^ = 
p{x+L)-p{x) _ ^Qj^g^ Using the Reynolds decomposition 
of velocity field v = Ui + u where v ~ U {y)i, the Navier- 
Stokes equations for incompressible fluid (p = const, V • 
u = V • v) can be written as: 

d,iu,u) + Ud^u + UydyUi = - ^ + iiyd^U + uV'^u (2) 

and, since all derivatives dx"^ = dz'^ = 0, where is the 
mean value of an arbitrary flow property \1/ and z is a 
coordinate in the span-wise direction, we have: 

dxP/P - dyTa:,y = udyll (3) 

where the Reynolds stress —UxUy = Tx,y The relation (3) 

expresses the Reynolds stress in terms of an unknown ve- 
locity distribution U {y) . To close the problem, one has to 
write a differential equation for Tx^y, which involves a new 
unknown function, for which one has to derive another 
equation and so on ad infinitum. The procedure, leading 
to an infinite chain of partial differential equations, can 
easily be formally written down but is too hard to solve. 
Here we propose a global approach not relying upon in- 
formation about local properties of the flow. 
Integrating (3) in the interval < y < H yields 
Hdxp/p = —ul = vdU\H and denoting Txy = r and 
the centerline velocity Ucl = U{y = 0), we derive: 

Finally, introducing dimensionless parameters Z = 
y/H = y+/R*, y+ = yu^/v, U+ = U/u*, i?* = u*H/u 
and Tx^y = M*T+, the equation valid for both laminar and 
turbulent flows is readily derived: 



U-UcL = - 




= ~R,-ua! dZ [ T+{Z')dZ' (4) 
6 Jo Jo 



To evaluate the integral (4) we need an expression for 
T+(y), which, at this time, is impossible to derive with- 
out generating an infinite chain of partial differential 
eqautions. Instead, let us define a thickncs of sublayer 
y = VsL, which combined with the expression (4), gives 
an exact magnitude of a global property U — Ucl- In 
other words, the integral in the right side of (4) is: 

U-UcL = -^R*+u,R, dZ Z'dZ'-I 

6 Jo Jo 

with 

/ — Ut^R^ / ^ dZ T+{Z')dZ' — Ut:y^]^ < t+ > 
Ji-^ Jo 

and < >< denoting the mean value of dimensionless 
Reynolds stress in the sublayer R^ — y'^j^ < y+ < R*. 
Thus, 

U-UcL = -au* + 0{1/R,) (6) 

where a = yfi^{^+ < t+ >)• The formula (6) defines the 

Taylor expansion in powers of a small parameter u^/U. 
We can see that as i?* — > cxd, the dimensionless parameter 
ih = 1 — oc — > which reflects the fact that, with 

/ UcL UcL ^ ' 

increase of the Reynolds number, the velocity profile U{y) 
flattens. This small parameter is crucial for the theory 
developed below. Based on numerical and experimental 
data y^j^ w 30, and | < t+ > | w 0.4 - 0.45 gives a « 4. 
(See, for example Ref. [4]),. 

The formally exact relation (6) has recently been verified 

by Zagarola ct al [5] in experiments on the Princeton 
SuperPipe giving a « 4.3 for 10^ < Reo < lO*". Sim- 
ilar result can be obtained by integrating the relation 
UcL - U{y) = u^F{^) [6] in the interval < y < R. 
This gives Ud — U = au^ where a = F{x)dx. The 
shape of the function F{x), consistent with logarithmic 
velocity distribution, was used by Prandtl who, based on 
experimental information, obtained a ~ 3.75. The later, 
probably more accurate, measurements gave a ~ 4.0 (see 
Ref. [7] and references therein). 

It follows from (6), that the skin friction in the pipe flow 
is equal to: A = 8(f f = ^(^^ - l)^ 

The predictions from this relation with \ w 0.42, are 
very close to experimental data collected from a smooth 
pipe by McKeon et al and from honed and commercial 
rough pipes studied by Schockling et al and Langelansvik 
et al [8], respectively, which is close a w 4. estimated 
above. 

As i?e ^ oo, the sublayer dominated by intermittent 

bursts of velocity derivatives, dissipation and produc- 
tion, can be considered as a low - Re turbulent flow with 
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the mean velocity U{ysL) oc and the y-component of 
the fluctuating velocity w{ysL) oc w*. The kinetic en- 
ergy generated in the subalyer (0 < ?/+ < 15 — 30) is 
rapidly mixed and dissipated in the bulk. It is clear that 
in the large Reynolds number limit ysL — ^ 0, the mean 
energy flux through the separating surface {x,ysL,z) is 
pK{ysL)w{ysL)LW and the energy balance can be writ- 
ten as : 



pK{ysL)w{y,L)LW = 0{uILW) = pSLWH (7) 

where £ = ^ j" £{y)dy = 0{ul/H). 
Flat Plate Boundary Layer. We consider a flat plate < 
a: < oo and y = 0- The freestream velocity of incoming 
flow is Uo = Uq'i and we are to analyze the Navier-Stokes 
-Prandtl equations in the boundary layer approximation: 



dU dV 
dx dy 



= 



U— V— - —iv— T 
dx dy dy dy ''^ 



(8) 



(9) 



As a; — >■ oo we, assuming self-similarity of the velocity 
proflle write: U = U{^) = U{r,), V = F(^) = V{v) 
and Tx,y = Tx.y{j^) = Tx,y{ri) where the defined below 
width of the boundary layer 5{x) must be found from 
equations of motion. The incompressibility constraint (8) 
gives: 



V{x,y) = - f 

Jo 



dU{x,y'] 
dx 



d6{x) r\,dU{r]')^ 

'0 



dx 



drj' 



dy' = T / v'—nr^dv' 

(10) 

Integrating (9) over the interval < y < oo, and in- 
troducing the 'displacement thickness' 6 we, using (10), 
express the skin friction in terms of the boundary layer 
thickness 5: 



dx 



d5{x) 
dx 



f 

Jo 



Uo Uo Uq 



where u1 = z^^lo- Based on as yet unknown function 
S{x), we define an averaged- over- the -boundary -layer 

property * * ~ j Jo"* ^{y)dy = ^{ri)dri. Since at the 
edge of a boundary layer y = d{x), the velocity is J7 = 
U{d{x)) and kinetic energy K = K{d{x)), the familiar 
integral balance equations must be somewhat modified. 
For example, integrating the differential energy balance 
equation: 



in the interval < y < 6{x) and recalling that 
dK dV 
we, using an incompressibility constraint, derive: 



f 

Jo 



dx 



K{x, y)U{x, y)dy + V{x, d)K{6) = 



p6 PiTT 

■ / T^y—dy-5£ + Q{5) (12) 
Jo oy 



where Q{5(x)) = w{5)u^(5) is small. In the limit ^ ^ oo, 
V{x, 5) — >■ and the relation (12) tends to a famil- 
iar energy balance (sec for example Hinze [xx]). With 
V{x,5{x)) = -£((5C7) + ^U{5), directly following from 
(10), the relation (12) takes a very simple form: 



dU - 
Txy-Q^dy - 6£ + Q{5) 



(13) 

Now, based on the results of a previous section, we in- 
troduce a new definition of a boundary layer thickness: 



K « bul: 



K{S) 



(14) 



Since Uq > U{5) > U, then, according to proposition 
(W, l-^ = i- ^+«^>0. Therefore, 



U{5) — U = 0{u^,), which crucially differs from a widely 
accepted ad hoc engineering definition of the boundary 
layer thickness U{5) = 0.99C/o implying tp = 1 — U /Uo = 
0(1). Let us demonstrate that the anzatz (14), combined 
with the energy balance (13), leads the well-known em- 
pirical relation A oc j^^ttj^t^- As follows from (14) and 
(11), 7fC7 w /3m2(J7o - au.,) and: 



dx 



uiUn 



0{fj « ul; 



dS_ 
dx 



dx 





(15) 



It will become clear below that as x — oo, 
0{uI/Uq) . Substituting the anzatz (14) into the energy 
balance (13) and accounting for the estimates (15) we, 
equating the terms of the same powers in and neglect- 
ing the O{u^/Uo) contributions, derive (3 = c, 



f 

Jo 



dU 



dy = S£ + Q{S) = 0{ul 



and 



UoO-T— = aiu^ — 
dx dx 



(16) 
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where ai = a — a < is an unknown parameter of this 
global approach which must be obtained from a full local 
theory. It will become clear below that the balance (13) 
is possible only if ai < 0. 
It is easy to see that the expression: 



In^^' 



, dd 
A oc — 
ax 



(17) 



with K = S/af is a solution to (16). Indeed, integrating 
(16) and dividing the outcome by Uq, we obtain: 



A 
2 



lailv^ f dX 



4\/2 



/ 



dx 



dx = 



4\/2 



IS a 



This result shows that the anzatz (14) with A oc 
solution to the Navier-Stokes -Prandtl equations of mo- 
tion. 

Setting for a time being all proportionality coefficients 
equal to unity, we introducing Sq = Res = and 

Rsx = solve the the differential equations (17) with 
the result: Res[{ln Reg)'^ — 2hi = Re^ and, as Res — >■ 

DO, 



S{x) 



hi' 



(18) 



So 



Summary and discussion. 1. In this paper, based on the 
Navier-Stokes equations, for a channel/pipe flow we de- 
rived the scaling relation (6), valid in pipe/channel flows 
the large-Reynolds number limit. Known for many years, 
this formula was previously obtained from analysis of 
experimental data or assumed scaling of velocity profile 
U{y). While the linear scaling with is an exact con- 
sequence of the Navier-Stokes equations, the magnitude 
of parameter a w 4. was estimated in this work from the 
near - wall data on the Reynolds stress. 

2. This result led to a new dynamic definition of the 
boundary layer thickness given by (14). 

3. The expression (17) for friction factor (skin friction) 
was found as a solution to the Navier-Stokes-Prandtl 
equations. 

4. These results are accurate up to the O(u,/J7o)- 
corrections. 
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